Effects of non-uniform acceptance in anisotropic flow measurement 
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Applicability of anisotropic flow measurement techniques and their extension for detectors with 
non-uniform azimuthal acceptance are discussed. Considering anisotropic flow measurement with 
two and three (mixed harmonic) azimuthal correlations we introduce a set of observables based on 
the X and y components of the event flow vector. These observables provide independent measures 
of anisotropic flow, and can be used to test self-consistency of the analysis. Based on these observ- 
ables we propose a technique that explicitly takes into account the effects of non-uniform detector 
acceptance. Within this approach the acceptance corrections, as well as parameters which define the 
method applicability, can be determined directly from experimental data. For practical purposes a 
brief summary of the method is provided at the end. 

PACS numbers: 25.75.Ld, 25.75.-q, 25.75.Gz 



I. INTRODUCTION 

One of the most important observable in the study 
of ultra-relativistic nucleus-nucleus collisions is the 
anisotropic transverse flow P, [3, H, [l]- It is quanti- 
fied by coefficients in Fourier decomposition of particle 
azimuthal distributions relative to the collision reaction 
plane 0, Q , the latter is defined by the beam direction 
and the impact parameter. The second harmonic in such 
a decomposition is called the elliptic flow. Observation of 
strong in-plane elliptic flow increasing with collision en- 
ergy from top AGS (Alternating Gradient Synchrotron) 
energies Hp, then at CERN SPS (Super Proton Syn- 
chrotron) [3, and recently at RHIC ('Relativistic 
Heavy Ion Collider) [ul, [H, M, M, M, M, [H shows 
that the bulk matter created in the high energy heavy 
ion collision strongly interacts and behaves collectively. 
Taken together with a key feature of anisotropic flow to 
be sensitive to the early stage of the collision, this indi- 
cates rapid attainment of local thermal equilibrium in the 
created system. For central collisions, RHIC results are 
in agreement with ideal (zero viscosity) hydrodynamic 
predictions [H, [l^, H^l . These observations have played 
an important role in the discovery of the strongly inter- 
acting Quark Gluon Plasma (sQGP), the new form of 
matter formed in heavy ion collisions. 

There exist a set of different techniques for anisotropic 
transverse flow measurement, which have been success- 
fully applied at a variet y o f experimental setups world- 
wide 0, a i, M, M MM, MM [H, [13 • As the reac- 
tion plane is not known experimentally, various methods 
exploit azimuthal correlations between particles as all of 
the particles are correlated to the same reaction plane. 
The most widely used are the event plane m M, t he scalar 
product [2l|, [l^], and mixed harmonic [a. l23l. [23 . [25j 
methods. Analysis of anisotropic flow with azimuthal 
correlations requires to examine contribution from effects 
not related to orientation of the reaction plane, such as 
resonance decays, jets, as well as effects of conservation 
laws (charge, momentum, etc). Most of these so called 
non-flow correlations are due to few particle interactions. 



and their relative contribution scales inversely propor- 
tional to particle multiplicity in the event. For that rea- 
son in multi-particle correlations these effects are sup- 
pressed compared to collective effects such as anisotropic 
flow. The complete consideration of methods to study 
non-flow effects is beyond the scope of this paper (for 
more discussions see, for example, [2^ and references 
therein) . In this paper we assume that the azimuthal dis- 
tributions of particles produced in the collision depends 
only on the orientation relative to the reaction plane. 

High statistics experimental data collected in recent 
years at RHIC allow to p erform very precise measure- 
ments 0, i i, M, M, M M M MmE^- with the 

new data, the systematic uncertainty in the analysis be- 
comes dominant compared to the statistical errors. It is 
vital to carefully investigate the systematic errors, in par- 
ticular due to non-perfect azimuthal acceptance, as well 
as to review the applicability of different methods in this 
respect. In addition, the availability of different exper- 
imental setups with complicated azimuthal acceptance, 
such as central arms of PHENIX (A Physics Experi- 
ment at RHIC) detector [MUS, NA49 TPC (Time Pro- 
jection Chamber) [H, and PHENIX/STAR (Solenoidal 
Tracker At RHIC) ZDC SMDjZero Degree Calorimeter 
Shower Maximum Detector) [231, requires further devel- 
opment of new and/or generalization of known methods 
for use with detectors with significant acceptance non- 
uniformity. Such generalization would allow to enrich 
already available data with anisotropic flow measure- 
ment results from a wider range of experimental setups. 
This further provides an opportunity for a comprehensive 
comparison of available theoretical predictions against 
the experimental data. 

In this paper we describe a procedure to broaden the 
applicability of known methods to measure anisotropic 
flow for a range of detectors with non-uniform azimuthal 
acceptance. Distinctive feature of the suggested ap- 
proach is that the acceptance corrections can be deter- 
mined directly from experimental data. This eliminates 
the need to perform time-consuming and model depen- 
dent Monte-Carlo simulations of the detector acceptance 
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and efficiency. The main idea is demonstrated on an ex- 
ample of two-particle correlations, but for completeness 
we also provide formulae for the case of three particle 
correlations (mixed harmonic) technique, which, as dis- 
cussed above are less susceptible to non-flow effects. We 
keep the discussion of more complicated three particle 
correlation case in separate subsections, such that they 
can be easily skipped if the reader is interested only in 
the main idea. 

The paper is organized as follows. In section [IT] we re- 
view notations and definitions of basic quantities used in 
anisotropic flow analysis. We formulate them in a way 
that later helps us to introduce new observables based 
on different components of the event flow vector. These 
observables provide independent measures of anisotropic 
flow, and can be used to test self-consistency of the re- 
sults. In section Hm we discuss the effects of non-uniform 
detector acceptance, and describe the procedure of cor- 
recting the single particle and event flow vector such 
that observables derived for the perfect detector can be 
used. We provide the smallness parameters that can be 
used to quantify the range of applicability of the pro- 
posed method. These parameters can be estimated di- 
rectly from experimental data. Finally, we summarize 
the method. 

A similar problem of flow analysis with non-uniform 
acceptance detector was discussed in [s^ [3l| in the con- 
text of cumulant and Lee- Yang zeros analyses based on 
the use of generating functions. Though in some cases 
these techniques might yield to similar results, our inde- 
pendent approach clarifies the nature of the acceptance 
effects in flow studies, and further provides the required 
tools for the most often used analyses based on correla- 
tions with the event flow vector. 



II. METHODS 

A. Definition and notations 

Anisotropic transverse flow of particles produced in 
heavy ion collision is quantified by coefficients in Fourier 
decomposition of particle azimuthal distribution In 
this paper we use particle azimuthal spectra normalized 
to unity (particle production probability density): 



P{(t) - *Rp) 



\ n—l 



w„cos[n(0- ^'rp)] . (1) 



Here (j) is the particle azimuthal angle, and w„ is the n- 
th harmonic anisotropic flow. The reaction plane angle 
in equation ([1]) is an azimuthal angle of the impact 
parameter. In general, coefficients Vn are functions of 
particle transverse momentum pt and pseudorapidity rj. 
Such dependences are hereafter assumed implicitly and 
not indicated in the notation for the sake of brevity and 
simplicity. 



We introduce a unit vector ?i„ (n > 0), defined for each 
particle: 

Un = Xn + ii/n = COS + i siu 710 = exp{m(/)}. (2) 

In this equation ?i„ is given as a complex number with an- 
gle n(j) and absolute value of unity. Throughout this pa- 
per we use complex number representation, but all equa- 
tions can be re-written in terms of commonly used alge- 
bra of 2-dimensional vectors. An estimate of the reaction 
plane orientation is usually obtained with the so called 
ri-harmonic event flow vector Qn, which is defined as a 
sum of w„-vectors over a specific subset "EP" of particles 
produced in the collision: 



)„ = ^Un= ^( 



cos 710 



I Sm 71(7 



EP 



EP 



Xn + iYn^\Qn\ exp{m*5^p}. 



(3) 



Here X„ and Yn are the event vector components, and 
"I/gp is n-th harmonic event plane angle. For each n the 
event flow vector Qn (or 5'ep) provides an independent 
estimate of the reaction plane orientation. Note that Un 
itself can serve as a (5„-vcctor. However this is not very 
practical, since the more particles used to define the event 
plane the closer its orientation will be to that of the re- 
action plane. 

To derive the main formulae of proposed technique we 
find it very useful to introduce, though experimentally 
unobservable, u„ averaged over events with fixed orien- 
tation of the reaction plane: 



(^")*RP + * (yn)*RP = / - *RP) 



27r 



1 + 2 ^ cos [m {(j) - ^fpp) 



Vn (cos n^'Rp + i sin Ti^fRp ) . 



(4) 



In this section we consider the case of a detector with 
perfect azimuthal acceptance. This implies that an inte- 
gration over azimuthal angle (jj goes over 27r without any 
weight. We will relax this assumption in section|TTT]when 
considering detector acceptance effects. Similarly: 

= (Mvn) EP (cos 774'rp + i sin n^Ppp ) 

= T/„(cos7i\I'rp + 7sin775'Rp). (5) 

Here Vn = (Mu„)gp is an average 7i-harmonic anisotropic 
flow Vn convoluted with multiplicity M of particles from 
a subset "EP" used to calculate the event flow vector Qn- 
From equation ^ it follows that Qn is defined as a 
vector in transverse plane, which on average has an ori- 
entation of that of the reaction plane. This feature can 
be used to define the event flow vector with detectors 
without tracking that are only sensitive to the shape of 
the particle distribution in the transverse plane (for ex- 
ample, calorimeters). The only, but very important, re- 
quirement to be fulflUed is that the Qn components, X„ 
and Yn, should be on average proportional to cos 77 5* rr 
and sin 71 5* Rp, respectively. 
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B. Anisotropic flow from diff'erent components 

1. Two particle correlations 

Anisotropic flow via two particle correlations can be 
obtained with the so-called scalar product method [2^ . 
According to the scalar product technique, one considers 
the average of the product of u„ and Q„ vectors over all 
events. With the help of Eqs. (jl]) and ([SJ this average 
can be written as an average over all events with fixed 
reaction plane with further average over all reaction plane 
orientations: 



d^ 



R,P 



27r 



{UnYn 



(6) 



to (cos [n{(f) — ^'ep)]): and Eq. ([8]) leads to the main ob- 
servable of the event plane method @: 



cos 



cos 



"(*EP - *EP) 



(12) 



Similarly, the second equality in formula (1111) gives an 
observable used by the NA49 Collaboration |10l |: 



^/2 



sinn0 • smnWgp 



smnWgp • smnW^p 



(13) 



Here and in Eq. wc use the event plane angle defined 
for the subevents, which resulted in an extra factor of two 
compared to Eqs. ((51 [TT|) . 



Here, angle brackets with subscripts 4'rp, {...)^^^, de- 
note the average over events with fixed ^'rp ; angle brack- 
ets without subscripts, (...), correspond to the average 
over entire ensemble of events with all possible orienta- 
tions of the reaction plane. 

The left hand side of Eq. ^ can be measured from 
experimental data. To obtain w„ one needs to evaluate 
Vn- This can be done by using random sub-events, i.e. 
randomly assigning particles used to construct the event 
flow vector into two subsets a and b 0: 

(QIQ':) = + {Y^Y^) = \vl (7) 

The factor of 1/4 here takes into account the multiplicity 
difference between the full event and that of sub-events 
a and b. 

From Eqs. ^ and © we obtain: 



(8) 



We further note that the two terms in Eqs. (O and (O are 
statistically independent, which allows to consider them 
separately: 



{XnXn) = {VnYn) = ^"tiK, 



(9) 
(10) 



thus providing two independent measures of anisotropic 
flow: 



{ynXn) 



(11) 



Independent observables (fTT|) can be used to check the 
self-consistency of the results. 

Note, that with normalization of the (5„-vector to 
unity, Q„ Qn/IQnl, the average {u„Q*„) in © reduces 



2. Three particle correlations 
In the case of three particle correlations one considers: 

= {xlxiX^n - v'^.vlX^n + xlv^^Y^n + V^iY^n) (14) 

= / ^ «)*HP «)*.P (Q2n)*,p - VIV^^, 



where: 

Then, the anisotropic flow u„ is given by: 



(15) 



(16) 



If one normalizes (5„-vector to unity, formula (1161) l eads 
to an observable of the mixed harmonic method [21, : 



<^cos 




I) 




2n(*EP - *EP ) 


) 



(17) 



All terms in formula ([14)) are statistically independent, 
which leads to a set of equalities: 

= HV2n. (18) 

Thus one obtains four independent observables to mea- 
sure anisotropic flow from three particle correlations: 



V2 



{xlxiX2n) 



-V2- 



{v^y';,X2n) 



V2 



{KyiY2r. 



^/2^^M^. (19) 
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As in the case of the two particle correlations, each of 
the four terms in Equation provides an independent 
measure of anisotropic flow, and can be used to check the 
self-consistency of the results. 

Note, that one can construct three particle correla- 
tion function from u„ and Qm vector components other 
than that defined by formula p4p . Some examples are 
{unQnQ2*n) {QnQn'^2n) ■ Derivation of observables 
based on these correlators are similar, but in this pa- 
per we only consider combination (jl4p , which in the case 
of rt = 1 leads to the known observable for directed 

flow [13, il. 



III. EFFECTS OF NON-UNIFORM 
ACCEPTANCE 

In order to generalize our consideration for the case 
of imperfect acceptance we introduce the acceptance 
function ^(0) which we normalize to unity (similar to 



(20) 



Then the average of some function /(0), which depends 
on particle azimuthal angle (j), at fixed reaction plane 
orientation is given by the integral: 



(/)*HP = j dM{<f>)fi^)p{^ - *RP) (21) 

oc 

= f + 2^v,n [/c„, cos m*Rp + /s„ sin m*Rp] . 

m— 1 

Here for brevity we introduce notation c,„ — cos m(j) and 
Sm = sinmt/), and denote by /, the average over the 
detector acceptance: 



2tt 



A{cb)f{cl>). 



(22) 



One might note that c„ and s„ represent n-th harmonic 
coefficients in the Fourier expansion of the acceptance 
function A((f>). An important observation is that the ac- 
ceptance average of /, /, coincides with the event aver- 
age, (/): 



(/) = 



J d'ifRpd<j)A{(f>)p{(f> - *Rp) 



/■ (23) 



We assume here that the distribution of the reaction 
plane angle, ^'rp, is uniform within a given centrality 
event sample. Experimentally, this can be achieved by 
using for the collision centrality determination the inde- 
pendent detector with a good azimuthal coverage. Con- 
sequently, all acceptance average quantities can be ex- 
tracted directly from the experimental data by the corre- 
sponding average over all particles in the event sample. 

Formula (|2ip allows to re-write the expressions for 
(xn)^ and {yn)^ taking into account the effects of 



the non-uniform detector acceptance. For clarity of com- 
parison with Eq. ^ we separate the term with m ^ n, 
which is the only non-vanishing term in case of perfect 
acceptance: 

(x„)^^p = c„ -I- u„a^„ {cosn^Rp + sinn*Rp 

00 

+ J2 ([^"-™ + C0STO*RP 

(j/n)^j,p = s„ + w„a^„ {sinn*Rp + Aj" cosn^-Rp 

+ {l^n-m - K+ni] siuTO^Rp 



+ Ar'-m]cosm*Rp)}.(25) 



Here we have introduced the acceptance coefficient 



afn = 1 ± C2„ = 1 ± cos 2n(/), 



and the following smallncss parameters: 



„,^.„ — T — : A„,^ 



(26) 



(27) 



These parameters define the relative contribution of dif- 
ferent terms in Eq. ([24|) and ([25|) . For a particular case 



of 771 = n, values of Aj^^ are defined only by detec- 
tor acceptance, while in general they also depend on the 
ratio of anisotropic flow Vm and u„. For a perfect de- 
tector Cm = Sm. = (and consequently all parameters 



c,sib 



0), and Eqs. ^ and ^ are reduced to Eq. (g]). 



Eqs. ([24]) and ([25|) show that acceptance effects result 
in coupling of equations for flow of different harmonics, 
and in general a simultaneous analysis of all harmonics 
is required. However we proceed below neglecting m ^ n 
terms. The relative contribution of m ^ n terms is of 
the order of Aj^'^^., and this case can be understood as 
either when the n-th harmonic flow is dominant: w„ 3> 
Vm^n (such an assumption, for example, is made in [31 1 
when discussing acceptance effects), or the acceptance 
effects for higher harmonics are small: c, s„zp„i/a^,j <C 1, 
or both. 

In the following, we distinguish three types of accep- 
tance effects: 

1. Shift of the Un-vector due to non-zero values of Cn 
and s„ in Eqs. (pi)) and (pS)) . This effect can be cor- 
rected for by subtracting from the w„-vector com- 
ponents their corresponding average values: 



Xn ; Un Vn -^n ■ 



This procedure is called re- centering. 



(28) 
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2. Twist of the Un-vector that results in appearance of 
sinn\l/Rp (cosnVpRp) terms in Xn [yn] components 
of it„-vector. Determined by non-zero values of 
in Eqs. and (PS)) . this effect can be corrected for 
by the diagonalization procedure (after re-centering 
has been applied): 



1 



Vn 



y 



^2n 



1 



^2n ^2n 



(29) 



Twist effect is zero if sin 2n(f) ~ 0, e.g.. in case n = 1 
for detectors symmetric in x or y (such as of rect- 
angular shape). 

3. Rescaling of the u„-vector, which is defined by the 
coefficients a^„ in Eqs. ([24]) and {251). This effect 
is the most important one next to the shift of the 
M„-vector, and it can be corrected for by rescaling 
the u„-vector components with acceptance coeffi- 
cients (after the re-centering and twist correc- 
tions have been applied): 



^2n 



Vn 



'■2n 



(30) 



The acceptance corrected u"'-vector has the same aver- 
age (it'")^^p as in the case of a detector with perfect 
acceptance (compare with Eq. 



= u„(cosnvI'r{p + i sinn\I>Rp). (31) 



Similar corrections can be applied for the Qn-vector com- 
ponents, which we write in the following form (contribu- 
tions from m ^ n terms have been neglected): 

(^«>*aP = + (cosTi^-Rp + sinn^Rp) , 
(^»)*RP = F„ + A2"„(sin7i*RP + A^;^cosn4'Rp). 



^2„A2„- Applying 



The symmetry requires ^J„A2^ 

corrections (|28ll30p for the Qn-vector (a^„, have to 
be replaced with A^„), one gets: 



(<?"') *Rp = C0S7ivI'rp + isinn^'Rp. 



(32) 



From Eqs. (|3ip and (|32p it follows that all equations given 
in section Hi B I can be applied to u"' and Q"' vectors, and 
the same observables (|TT|) , (fT8|) , and ((T9|) can be used for 
acceptance corrected wj" and QJ" vector components. 

To clarify better how the above described corrections 
work we consider below the correlations between uncor- 
rected u„ and Q„ vector components, and discuss what 
kind of effects are removed by each of the corrections 
(12811301). 



A. Two particle correlations 

Acceptance effects in conjunction with anisotropic flow 
may lead to various spurious correlations, such as corre- 
lations in multiplicity and/or transverse momentum. In 



particular, multiplicity correlations in two kinematic re- 
gions a and b arc given by the following equation: 



[PaPb] = J 



RP 



2tt 



d4)ad(l)bA{4>a)A{4>b)paPb 



1 + 2 



oo 

E 

rn—1 



\{v )J, (33) 



where u,„ = Cm + iSm, and denotes the av- 

erage over kinematic regions a and b. According to 
Eqs. [55]) . the (event) average of the product of two 
functions / and g defined in regions a and b can be writ- 
ten as: 



ifagt) 



1 



[PaPb] 



RP 



2tt 



(fa) 

*RP (9b}<s,„ 



(34) 



Deviation of the denominator from unity (which is the 
value for no multiplicity correlation) is defined by non- 
zero terms in the sum in Eq. p3p over m-harmonics. 
Taking into account that the measured anisotropic flow 
at RHIC is Vm < 10%, in the most pessimistic estimate, 
using {um)a fc Ij we obtain |[pa/06] — 1| < 2%. In prac- 
tice, {um)a 6 < 0-2 + 0.3, which reduces the acceptance 
effects on multiplicity correlations to the level of a tenth 
of a percent. In principle, such effects can be consistently 
taken into account, but for the sake of simplicity, below 
we proceed neglecting these multiplicity correlations. 

In analogy to ([6]) we consider the correlations between 
uncorrected m„ and Qn vector components: 

(XnXn) ~ XnXn + -^a2n^2ri + ^2n^2ri) ' 

{VuYn) = y„F„ + ya2,,A2-„ (1 + A^;A^;) . (35) 

The first terms, c„X„ and SnYn, can be removed by 
the re-centering procedure ([28)) . Note, that c„ and X„ 
enter as a product, what allows to re-center only the 
event vector components. Similar, to remove the second 
order terms, A2^A2^. it is sufficient to apply the twist 
correction only for the Qn-vector components. Param- 
eters and A2^ can be obtained with the random 
subevent technique. In that case, they are defined by a 
set of coupled equations: 



(36) 



^2n 



After the re-centcring and twist corrections have been 
applied Eqs. ([55)) leads to the following observable for 
the anisotropic flow «„: 



1 



1 



al, ^2 {X--X'^) a2„ VWFW) 



(37) 



Rescaling of the M„-vector reduces this to Eq. ([TT)). 
which should be written for rescalcd u„ and shifted 
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and twisted Q" vectors. Note, that correction factors 
yJ{X'^°-X'^'') and ^/Wll^^j!^ correspond to rescahng of 
the Q"-vector. 

Eq. ((57)) shows that in case of two particle correla- 
tions re-centering and twist corrections of the M„-vector 
are not required. This equation can be also used for the 
QJ^-vector, if the second order corrections defined by the 
terms Aj^Aj^ are small, and can be neglected (twist cor- 
rection is not required). 

B. Three particle correlations 

In the case of three particle correlations we consider: 

{xlxiX2n) = «4) X2n + K (4^2„) + 4 {<X2n) 

+ ^-f4n'At (l - Kf + 2A^+A1+) (38) 

Similar expressions can be written for other terms in 
Eq. (dH]). In contrast to the case of two particle correla- 
tions, re-centering procedure (|28p is required for all three 
vectors u^, u\, and Q2n- Acceptance coefficients and 
A4^ are given by a set of coupled equations p6|) written 
for the (52n-vector. Twist corrections applied to Q2n and 
vectors removes the terms X?^ and 2A2^A4^, and 
Eq. ([38]) leads to observable ([HI) written for u'^ and Q"' 
vector components. 

IV. METHOD SUMMARY AND CONCLUSION 



Non-uniformity of the detector acceptance is quantified 
with coefficients c„ = cosncj) and s„ = sinn0, which fur- 
ther define coefficients and A^!^ given by Eqs. ([26]) 
and ([?7|) . Though accounting for the acceptance effects in 
general might be difhcult as it requires a solution of a set 
of coupled equations with different harmonics involved, 
in the case when the contribution of m 7^ n harmonics 
can be neglected, the problem significantly simplifies. It 
becomes possible to correct the single particle m„ and 
event flow Qn vectors such that the conventional observ- 
ables (derived for the perfect detector) can be used. Note, 
that both, the acceptance coefficients c„ and s„, and the 
correlators between it„ and Qn vectors can be obtained 
during a single pass over the data. This can significantly 
reduce the amount of time needed for the calculation. At 
the same time, due to variation of detector acceptance 
in time, with collision ccntrality, vertex position, etc., it 
may be important to apply the acceptance corrections 
separately run-by-run, for different vertex position, etc. 
In that case it might be more convenient to split the pro- 
cedure into a few steps with two passes over the data. 
During the first pass the acceptance coefficients c„ and 
s„ are extracted as a function of different ccntrality, time, 
etc. and all coefficients needed for acceptance corrections 
presented in Eqs. (|28ll30p . both, for m„ and Q„ vectors, 
are calculated. During the second pass over the data the 
correlators of the standard procedure given by Eq. ^ 
[To)) (Eq. in the case of three particle correlations) 
are calculated. Finally, the flow values are extracted as 
given in Eqs. HD), (flj)) . 



In this paper we discuss new observables for anisotropic 
flow measurement based on correlations of x and y com- 
ponents of the flow vectors. Providing independent mea- 
sures of anisotropic flow they can be used to check self- 
consistency of the analysis. Moreover, these observables 
allow direct accounting for acceptance effects, which we 
discuss in detail for two particular cases of anisotropic 
flow measurement with two and three (mixed harmonic) 
particle correlations. Importantly, acceptance correc- 
tions and parameters, which define applicability of these 
observables, can be determined directly from experimen- 
tal data. 
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